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ABSTRACT 



Most exoplanets have been discovered via radial velocity studies, which are in- 
herently insensitive to orbital inclination. Interferometric observations will show 
evidence of a stellar companion if it sufficiently bright, regardless of the inclina- 
tion. Using the CHARA Array, we observed 22 exoplanet host stars to search for 
stellar companions in low-inclination orbits that may be masquerading as plan- 
etary systems. While no definitive stellar companions were discovered, it was 
possible to rule out certain secondary spectral types for each exoplanet system 
observed by studying the errors in the diameter fit to calibrated visibilities and 
by searching for separated fringe packets. 

Subject headings: binaries: general — infrared: stars — planetary systems - 
techniques: interferometric 



1. Introduction 

In studies of exoplanet systems, certain assumptions are made about the inclination (i) 
of the systems discovered; i.e., it is assumed the orbit has an intermediate to high inclination 
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(i ~ 45 — 90°) because the probability of the orbit being nearly face-on (i ~ 0°) is extremely 
low. If we assume a random sample of orbital orientations along our line of sight, the 
probability of the inclination being less or equal to a given % is [1— cos(i)] while the probability 
of an inclination being greater than a given i is cos(z). Therefore the probability of an orbit 
with an inclination below 45° is ~30% while the probability of the orbit having an inclination 
above 45° is ~70%. 

This conclusion implies that the calculated companion mass, known only as the quantity 
m p sin 3 z where m p is the mass of the planet, is planetary in nature instead of stellar. While 
this is probably a safe conclusion for the majority of the exoplanets discovered, the chance 
remains that in a large enough sample, a few of the candidate planetary systems may be 
face-on binary star systems instead. 

We make no assumptions about the inclination. If a second star is present and is not 
more than ~2.5 magnitudes fainter than the host star, the effects of the second star will be 
seen in the interferometric visibility curve^. See Figure [1] for an example of the difference 
between the visibility curves for a single star and binary system. 

Two studies have shown that radial velocity observations of exoplanet systems alone are 
insufficient to distinguish between intermediate- to high-inc lination planetary systems and 



low-inclination binary star systems. IStepinski &: Black! (120011 ) estimated probability densities 



of orbital periods and eccentricities for a sample of exoplanet candidates and a sample of 
spectroscopic binary star systems with solar-type primary stars in order to determine if 
there were any fundamental differences between the two types of systems. They found the 
respective distributions of the two populations were statistically indistinguishable from each 
other in the context of orbital elements. 



In an earlier study, llmbert fc Prevotl (119981 ) modeled nine known exoplanet systems as 



binary star systems to test if the radial velocity observations could be explained by low-mass 
stellar companions. Although the probability of binary star systems appearing as planetary 
systems was low, ranging from 0.01 to 4%, the model results described the observations 
satisfactorily and showed it is possible for a binary star system to mimic an exoplanet 
system. 



Wu et al.l (120071 ) argue that ~ 2.5% of the exoplanet systems with Jupiter-mass planets 
with semimajor axes of <0.1 AU may have formed in binary systems where the binary star 
companion is highly inclined to the planet's orbit. Furthermore, it has been proposed that 



1 A limiting AK of 2.5 is a lower limit, as the true AK also depends on the absolute brightness of the 
two stars and could be higher for other systems. 
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planets form in the same manner in both single-star and binary-star systems during certain 
stages (ITsukamoto Making 120071 ). This would be significant due to the fact that more 
than half of all stars form in binary or multiple star systems. 



2. Interferometric Observations 

Our target list was derived from the general exoplanet list using declination limits and 
magnitude constraints. The stars needed to be north of -10°, brighter than V = +10 in 
order for the tip/tilt subsystem to lock onto the star, and brighter than K = +6.5 so fringes 
were easily visible. This reduced the original list to ~80 targets, and we observed 22 systems 
from 2004 January to 2007 September. 

The stars were observed using the Center for High Angular Resolution Astronomy 
(CHARA) Array, a six-element Y-shap ed interferometric array located on Mount Wilson, 



California (jten Brummelaar et al.ll2005l ). The Array presently uses visible wavelengths at 
470-800 nm for tracking and tip/tilt corrections and near infrared bands, H at 1.67 /zm and 
K' at 2.15 fim, for fringe detection and data collection. All observations of the host stars 
were obtained using the pupil-plane "CHARA Classic" beam combiner in the i^'-band. 

The observations were taken using many of the baselines the CHARA Array offers from 
intermediate-length baselines at 108 m to the longest baseline available at 331 m. Table 
[1] lists the exoplanet host stars observed, their calibrators, the baselines used, the dates of 
the observations, and the number of observations obtained. Calibrators were chosen to have 
small predicted angular diameters in order to reduce errors in the target's angular diameter 
measurement, and the calibrators were as close to the target stars as possible, usually within 
5°. This allowed for less time between calibrator and target observations, thus reducing the 
effects of changing seeing conditions as much as possible. 

We observed using the standard calibrator-target-calibrator pattern, which allowed us 
to see if the target star was changing with respect to the calibrator over time, assuming 
the calibrator did not have an unseen stellar companion or a circumstellar disk. In order to 
select reliable calibrators, we searched the literature for indications of binarity and performed 
spectral energy distribution fits to published U BVRIJHK photometry to see if there was 
any excess flux that would indicate a stellar companion. Calibrator candidates with variable 
radial velocities were discarded even if their SEDs displayed no characteristics of duplicity. 

Stellar companions in low-inclination orbits could show themselves in two ways: in the 
residuals to the diameter fit to calibrated visibilities (sensitive to companions at separations 
of ~0.5 to 10 milliarcseconds) and as separated fringe packets (for companions in the sepa- 
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ration range of ~10 to 50 milliarcseconds) . The two methods are described in more detail 
below. 



3. Characterizing Residuals to Diameter Fit 

The systematics in the residuals of the diameter fit to measured visibilities indicate 
whether or not a stellar companion may be present. For a single star, the residuals show a 
Gaussian distribution about (see Figure [2]). On the other hand, Figure [3] shows the strongly 
systematic behavior in the residuals for a known binary star system. Inspection of Figure [2] 
shows that the error estimates of the calibrat ed visibilities in all cas es exceed the (observed 



calculated) residuals (cr res ). As described in iMcAlister et al.l (120051 ). the error estimates in 



the measured instrumental visibilities are conservatively taken from the standard deviation 
from the mean of subgroupings of individual visibility scans. In this case, those estimates 
provide a reduced x 2 l ess than unity, which indicates that our individual visibility error 
estimates are too large. 

The standard deviation for the residuals for Figure [2] is ±0.056 while a rcs = ±0.106 
for Figure [31 The typical <7 res for a single star is well under ±0.100 while a res for a binary 
with a low AK is usually > ±0.100. F or example, the binary sy stem HD 146361, which is 



composed of two nearly-identical stars (jStrassmeier &: Ricdl2003l ) and has a AK ~ 0.2, was 



observed by Deepak Raghavan for three nights in May of 2007 along with the calibrator star 
HD 152598. The measured cr res for the three nights observations are ±0.180, ±0.229, and 
±0.118 indicating a departure from the single-star model for this system (private communi- 
cation). 

Another diagnostic to distinguish between the single star seen in Figure 2 and the binary 
star system shown in Figure 3 is to compare their x 2 values. A low x 2 value indicates a good 
fit to the single-star model while a large value indicates a poor fit, revealing the presence of 
a secondary star, circumstellar disk, or an asymmetry of the star. For example, the x 2 f° r 
the star shown in Figure 2 is 1.80 while x 2 = 20.89 for the binary star seen in Figure 3. 

For each exoplanet system, a variety of low-mass secondary stars were considered: G5 V, 
K0 V, K5 V, M0 V, and M5 V. Most of the stars in the sample are solar-type stars, so more 
massive main sequence spectral types need not be considered. The magnitude difference 
(AMk, listed as AK in the tables) and angular separation (a) of a face-on orbit between 
the host star and companion were calculated for each possible pairing: 



AM K = M s ± (m h - M h ) - m h , 



(1) 
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where Mh, s are the absolute magnitudes of the host star and potential secondary, respectively, 
nih is the apparent magnitude of the host star, and 

(m h - M h ) = 5 log (^j , (2) 

where 7r is the host star's parallax in milliarcseconds (mas). An estimate of the angular 
separation a in mas was calculated from Kepler's Third Law and %: 

a = [(M h + M s ) x P 2 } * x 7T, (3) 

where Mh jS are the masses (in M Q ) of the exoplanet's host star and potential secondary star, 
respectively, and P is the companion's orbital period in years. 

Additionally, the angular diameter 9 of the possible secondary was estimated using the 



calibration of radius as a function of spectral type from I Cox! ( 120001 ) and the parallax of the 



host star. The masses for the exoplanet host s tars a nd masses and radii for the possible 



secondary co mpanions were obt a ined f rom ICoxl (120001 ). which in turn were based on values 



derived from lHabets Heintzd (119811 ). who used observations of binary stars in order to 
create empirical relationships between various stellar parameters as a function of luminosity 
class. 

Tables El [3J and H] present the results of these calculations. For each exoplanet host star 
studied, the observed values required for the calculations described above are listed in Table 
[2J while Table [3] shows the calculated AK and a for each type of possible secondary star. 
Table H] includes the calculated angular diameters of the potential secondary companions 
for each system. It should be noted that if the companion star is a pre-main-sequence star, 
the resulting AK becomes smaller due to the star's increased brightness prior to hydrogen 
fusion and therefore has a higher probability of being detected. 

The resulting values for 9, AK, and a were then used to plot the visibility curves for both 
a single star with the host star's measured angular diameter and for a binary system with the 
calculated parameters. The projected position angle of a binary star vector separation onto 
the interferometric baseline is, of course, unknown and is here assumed to be 0° to explore 
the effects of the maximum separation exhibited by the secondary. As this angle approaches 
90°, the modulation to the visibility curve diminishes because the binary becomes unresolved 
at 90°. 

To estimate the detection sensitivity, the largest difference between the visibility curves 
for a single star and for a binary system with the parameters listed in Table [3] was calculated. 
This quantity, AV max , then represented the maximum deviation of the binary visibility curve 
from the single-star curve. 
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The lower limit to rule out stellar companions was selected to be 2cr res , where <7 res is 
the standard deviation of the residuals to the diameter fit; i.e., if AV max > 2cx res for a given 
secondary component, that particular spectral type can be eliminated as a possible stellar 
companion. If AV max < 2cr res , the effects of the companion would not be clearly seen in the 
visibility curve, and that spectral type cannot be ruled out. For each exoplanet host star, 
Table [5] lists the observed a res and the predicted AV mauX for each secondary type considered, 
and the final column indicates the cutoff point for the non-detection of a stellar companion. 
For example, if "K5 V" is listed in the last column, the spectral types more massive than 
a K5 V could be eliminated from consideration but stars of type K5 V and later are still 
possible companions. A dash in this column indicates all companion spectral types can be 
ruled out. 

Using the data from Table Figure H] was created to demonstrate the sensitivity of 
the interferometric observations to stellar companions. For each exoplanet host star ob- 
served, the absolute V- band magnitude was found from The Hipparcos and Tycho Catalogues 



(jPerryman et al.l 119971 ). Then the night with the lowest a ves was chosen to be the best case 



scenario when determining which secondary stars could be eliminated from consideration. 

The difficulty of ruling out the more massive companion types for the intrinsically 
brighter stars is expected, as the AK will already be large even for brighter companions and 
beyond the scope of the CHARA Array. Host stars that are less massive are fainter, and the 
AK lies more within the sensitivity limit of the CHARA Array. 

Two stars showed systematics in their visibility measurements that could indicate an 
unseen stellar companion in some datasets. The first star is v Andromeda e (v And, HD 9826, 



F8 V). Its first planetary candidate w as announced bvlButler et al.l (119991 ) and two more plan 
ets were discovered two years later (IButler et al.lll999l ). An M4.5 V stellar companion was 



found accompanying v And at a distance of ~750 AU from the central star (ILowrance et al. 



20021 ). This star would not have affected our search for more close-in stellar companions, 
as the angular distance from the host star is 55" and is well out of the field of view of the 
CHARA Array. 

v And was part of two intensive observing campaigns using the CHARA Array in 2005 
August and 2007 September and the data cannot be fit with a simple limb-darkened disk for 
two of the nine nights of data. In one dataset (2005/08/04), the target's visibilities briefly 
became higher than the calibrator's visibilities, and in the other dataset (2005/08/10), the 
target's and calibrator's visibilities separated over time. These patterns indicate one of the 
stars is changing with respect to the other. Two different calibrators were used for these two 
nights of data, and other data obtained using the same calibrators show no systematics in 
the visibilities. Therefore, we do not claim a stellar companion to v And at this time. 
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The second star to show oddities in its visibility measurements was p Coronae Bo- 



Noves e 


t al. 


( 


1997 


) before a 1 


planet ( 


Gatewood et al. 


2001) 



used high-dispersion infrared spectroscopy to determine if they could detect any flux from 
an M dwarf companion, as it would lie within the sensitivity limits of their instrument. No 
such flux was detected and they concluded the companion was planetary in nature. 

The controversy surrounding this system made it an interesting target to observe using 
the CHARA Array. We observed p CrB for three nights using three different calibrators 
and found that the data for two of the four nights exhibited behavior inconsistent with a 
single star. The two nights in question are 2005/06/29 and 2005/07/03, and because those 
observations were taken using the same calibrator and the data taken using other calibrators 
fall in line with a single star, it is likely that it was the calibrator and not p CrB that was 
varying. There are currently plans to i mage both v And and p CrB using the Michigan 
Infrared Combiner (IMonnier et al.ll2004j ) on the CHARA Array, which may help to clarify 
the situation for both stars. 



4. Separated Fringe Packets 

The second method to check for unseen low-mass stellar companions is by searching for 
separated fringe packets (SFPs). When a star has a wide companion (~10 to 100 mas), two 
fringe packets - one from each star - may be observed if the baseline orientation is favorable; 
i.e., if the projected baseline angle is approximately parallel to the position angle of the 
binary and both fringe packets are within the data collection scan window. However, if 
the two stars have a small angular separation or the position angle is perpendicular to the 
projected baseline angle, the two fringe packets will overlay each other and appear as one 
fringe packet. 

For example, if a binary system with a AK=0 has a separation of more than ~10 mas 
and is in the optimal orientation as described above, SFPs may be visible, as can be seen 
in Figure [5j If the same system has a AK ~2.5 or is not in the optimal orientation, no 
secondary fringe would be observed, as is shown in Figure El The baseline used in the 
observations also plays a role in whether or not SFPs will be detected, as the separation of 
the fringes depends on the baseline length. Therefore, a system appearing as an SFP on the 
long baseline of Sl-El will not be an SFP on the shortest baseline, S1-S2. 

The detection of SFPs also depends partly on whether both fringe packets lie within 
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the scan window. The width of the scan window depends on the baseline, wavelength used, 
and the frequency of the observations. For an average observation in the X-band using a 
100-m baseline, the scan window will cover ~300 mas while at a 300-m baseline, the scan 
window width is ~100 mas. If the SFP is wider than the scan window width, data on the 
second fringe cannot be collected. 

For completeness, all the stars observed in this project were checked for SFPs, whether or 
not the calculated separation of the secondary star would indicate the possibility of separated 
fringes. Each of the ~200 scans in every dataset consists of a sampling of the zero-path- 
length delay space where fringes are located. For each individual scan, the strongest fringe 
was located and an envelope was fit to the fringe. The fringe envelope was obtained using 
a Hilbert tranform, achieved by taking the Fourier transform of the fringe scan, setting 
the negative freque ncies to zero, and taking the modulus of the inverse Fourier transform 



(Born fc WoUHl959h . 



Then the peak of the primary fringe envelope for each of the data scans was located 
and shifted so that the fringe envelopes for all the scans overlaid each other and the fringe 
amplitudes were added together. This "shift and add" approach made it possible to view 
multiple fringe envelopes at once to see if there was any indication of a SFP. The result was 
a plot of the weighted mean fringe envelope. Figures [7] and [8] shows an examples of the fringe 
envelopes for a single star and binary system, respectively. 

If there were three or fewer bracketed observations in the dataset, only one observation 
from each of the calibrator and object stars was inspected for SFPs. Otherwise the first 
and last observations in a night's dataset were inspected from SFPs for both the object and 
the calibrator in order to maximize changes in the stars' position angles with respect to the 
baseline over time. No SFP binaries were discovered for the exoplanet host stars or their 
calibrators. 



5. Conclusion 

In an effort to cull out any possible binary star systems in the exoplanet sample, we 
inspected the exoplanet host stars using the CHARA Array using two methods. The first 
involved characterizing the residuals of the diameter fit to calibrated visibilities. If the 
observed residuals were more than twice the predicted variations for a given binary system, 
the secondary spectral type in question could be ruled out. No stellar companions were 
detected but certain secondary spectral types were eliminated for each exoplanet host star. 

We also inspected the data for secondary fringe packets, which could be present if a 
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secondary star had the proper orientation to the baseline used for the observations. No 
secondary fringe packet binaries were found, further reducing the possibility that these exo- 
planet systems are binary star systems instead. 

Many thanks to P.J. Goldfinger and Chris Farrington for their invaluable assistance in 
obtaining the data used here. The CHARA Array is funded by the National Science Founda- 
tion through NSF grants AST-0307562 and AST-0606958 and by Georgia State University 
through the College of Arts and Sciences. This research has made use of the SIMBAD lit- 
erature database, operated at CDS, Strasbourg, France, and of NASA's Astrophysics Data 
System. This publication makes use of data products from the Two Micron All Sky Survey, 
which is a joint project of the University of Massachusetts and the Infrared Processing and 
Analysis Center /California Institute of Technology, funded by the National Aeronautics and 
Space Administration and the National Science Foundation. 
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Table 1. Interfero metric Observations. 
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HD 




("TIT 1 ) 
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6 
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v And 
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W1-W2 (108 m) 
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11 
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14 
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10 
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15 
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13 
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5 






8671 
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16 
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20 
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5 
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22521 
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2007/09/11 
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2007/09/14 


5 
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8 










2005/12/16 
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6 
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2006/05/13 


10 
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2006/05/20 
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W1-W2 (108 m) 


2006/05/14 
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2007/02/05 


10 










2007/03/26 
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2007/03/30 


8 
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8 
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E2-W2 (156 m) 


2005/06/29 


6 










2005/07/03 


7 






146025 


W1-W2 (108 m) 


2006/05/12 


5 


177830 
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W1-W2 (108 m) 


2006/08/08 
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Table 1 — Continued 



Target 
HD 


Other 
Name 


Calibrator 
HD 


Baseline 
(length) 


Date 
(UT) 


Number of 
Observations 








Sl-El (331 m) 


2006/08/13 


6 


186427 


16 Cyg B 


184960 


Sl-El (331 m) 


2006/08/13 
2007/09/12 


6 
6 


190228 




190470 


W1-W2 (108 m) 
E2-W2 (156 m) 
Sl-El (331 m) 


2005/08/19 
2005/07/01 
2006/08/14 


5 
6 
8 


190360 




189108 


W1-W2 (108 m) 
Sl-El (331 m) 


2005/08/11 
2006/08/11 


10 

9 


195019 




194012 


W1-W2 (108 m) 
Sl-El (331 m) 


2005/08/11 
2006/08/06 
2005/10/23 


5 
6 
10 


196885 




194012 


W1-W2 (108 m) 
E2-W2 (156 m) 
Sl-El (331 m) 


2006/08/07 
2005/10/29 
2006/08/14 


10 

5 
5 


217014 


51 Peg A 


218261 


W1-W2 (108 m) 
Sl-El (331 m) 


2005/08/12 
2006/08/12 


14 

7 



Note. — The three arms of the Array are denoted by their cardinal directions: "S" 
is south, "E" is east, and "W" is west. Each arm bears two telescopes, numbered "1" 
for the telescope farthest from the beam combining laboratory and "2" for the telescope 
closer to the lab. 



Table 2. Exoplanet Host Star Observed Parameters 



Observed Planetary Parameters Calculated Values 



K n M s tar P K\ a p Afpsini a star sini f(M sta r) 

HD (mag) (mas) M Q (d) e (m s" 1 ) (AU) (Af Jup ) Ref (m-M) (AU) (M ) 



3651 


4.0 


90.0 


0.79 


62.2 


0.63 


15.9 


0.28 


0.20 


F03 


0.2 


7.06E-05 


1.21E-11 


9826 


2.9 


74.3 


1.3 


4.6 


0.03 


73.0 


0.06 


0.71 


B99 


0.6 


3.10E-05 


1.86E-10 


11964 


4.5 


29.4 


1.12 


2110 


0.06 


9.0 


3.34 


0.61 


B06 


2.7 


1.74E-03 


1.58E-10 


16141 


5.3 


27.9 


1.01 


75.8 


0.28 


10.8 


0.35 


0.22 


MOO 


2.8 


7.23E-05 


8.75E-12 


19994 


3.7 


44.7 


1.34 


535.7 


0.30 


36.2 


1.42 


1.68 


M04 


1.7 


1.70E-03 


2.29E-09 


20367 


5.0 


36.9 


1.17 


469.5 


0.32 


29.0 


1.25 


1.17 


B06 


2.2 


1.19E-03 


1.01E-09 


23596 


5.9 


19.2 


1.3 


1565 


0.29 


124 


2.88 


8.1 


P03 


3.6 


1.71E-02 


2.71E-07 


38529 


4.2 


23.6 


1.4 


14.3 


0.27 


53.8 


0.13 


0.77 


F01 


3.1 


6.81E-05 


2.06E-10 


59686 


2.9 


10.8 


1.15* 


303 





10.0* 




6.5 


M03 


4.8 


2.79E-04 


3.14E-11 


75732 


4.0 


79.8 


0.95 


14.7 


0.02 


72.2 


0.12 


0.84 


M02 


0.5 


9.73E-05 


5.71E-10 


104985 


3.3 


9.8 


1.6 


198.2 


0.03 


161 


0.78 


6.3 


S03 


5.0 


2.93E-03 


8.56E-08 


117176 


3.5 


55.2 


0.92 


116.7 


0.40 


318 


0.43 


6.6 


MB96 


1.3 


3.13E-03 


2.99E-07 


120136 


3.5 


64.1 


1.2 


3.3 


0.02 


469 


0.05 


3.87 


B97 


1.0 


1.43E-04 


3.54E-08 


143761 


3.9 


57.4 


1.0 


39.6 


0.03 


67.4 


0.23 


1.1 


N97 


1.2 


2.46E-04 


1.26E-09 


177830 


4.8 


16.9 


1.15 


391.6 


0.41 


34 


1.10 


1.22 


VOO 


3.9 


1.12E-03 


1.21E-09 


186427 


4.7 


46.7 


1.0 


800.8 


0.63 


43.9 


1.6 


1.5 


C97 


1.7 


2.51E-03 


3.29E-09 


190228 


5.4 


16.1 


0.83 


1146 


0.50 


91 


2.02 


3.58 


P03 


4.0 


8.30E-03 


5.81E-08 


190360 


4.1 


62.9 


0.96 


3902 


0.48 


20 


4.8 


1.33 


N03 


1.0 


6.29E-03 


2.18E-09 


195019 


5.3 


26.8 


0.98 


18.3 


0.03 


275.3 




3.51 


F99 


2.9 


4.63E-04 


3.95E-08 


196885 


5.1 


30.3 


1.27 


386 


0.3 


lO.Ot 


1.12 


1.84 


EE 


2.6 


3.38E-04 


3.47E-11 


217014 


3.9 


65.1 


1.0 


4.2 


0.01 


55.9 


0.05 


0.45 


M97 


0.9 


2.17E-05 


7.65E-11 



Note. — tHD 59686 and HD 196885 had no K\ given in the literature, so a K\ of 10.0 m s 1 was assigned for the calculations. 
*HD 59686's mass was obtained from lCoxl l l2000h . 

All K magnitudes are from the 2MASS All-Sky Catalog of Point Sources llSkrutskie et alj|2006h and 7r values are from The Hipparcos Catalogue 
iPerrvman et al.|[l997h. 

References: B97: iButler et alj Jl997h; B9 9 : iButler et all il999l); B0 6: IButler et all (1200611: C97 : ICochran et al.l Jl997ft; EE: Ex trasolar 
Planets Encyclopaedia ; F99 : IFischer et all lll999|); F0 1: IFischer et alj feOOl); F0 3: IFischer et alj 11200311: V197: iMarcv et all lll997 ); MOO: 



iMarcy. Butler, fc VoeH j2000h ; M02: IMarcv et alj 
N97: iNoves et alHl997h ; N03: iNaef et alj ll2003h : 



M03: iMitchell et al 
Perrier et alj d200: 



M04: 



I 20031) 
S03: ISato et al 



Mayor et all " 20041): MB 9 6: IMarcv fc ButleJ lll996h : 
20031 '): V00: Ivogt et al.l i bOOCh . 



Table 3. Predicted Parameters for Secondary Stars of Various Spectral Types 



G5 V KO V K5 V MO V M5 V 



Host 


AK 




AK 




AK 




AK 




AK 


, a . 


HD 




[ VYl O O i 
I 11 lcLS 1 




( m l 




[ m ac 1 




t lllda J 


l md 6/ 


( 111. <1S \ 


3651 


0.27 


1.15 


0.13 


1.12 


0.73 


1.09 


1.38 


1.05 


2.36 


0.97 


9826 


1.25 


0.22 


1.65 


0.22 


2.25 


0.21 


2.90 


0.21 


3.88 


0.20 


11964 


1.66 


12.83 


2.06 


12.55 


2.66 


12.28 


3.31 


11.90 


4.29 


11.12 


16141 


0.98 


1.37 


1.38 


1.34 


1.98 


1.31 


2.63 


1.27 


3.61 


1.18 


19994 


1.55 


5.32 


1.95 


5.22 


2.55 


5.12 


3.20 


4.98 


4.18 


4.69 


20367 


0.66 


4.75 


1.06 


4.65 


1.66 


4.55 


2.31 


4.42 


3.29 


4.13 


23596 


1.20 


10.81 


1.60 


10.60 


2.20 


10.39 


2.85 


10.10 


3.83 


9.51 


38529 


2.44 


0.48 


2.84 


0.47 


3.44 


0.46 


4.09 


0.45 


5.07 


0.42 


59686 


5.43 


3.53 


5.83 


3.46 


6.43 


3.39 


7.08 


3.28 


8.06 


3.07 


75732 


0.01 


0.45 


0.39 


0.44 


0.99 


0.43 


1.64 


0.42 


2.62 


0.39 


104985 


5.24 


2.84 


5.64 


2.79 


6.24 


2.75 


6.89 


2.68 


7.87 


2.55 


117176 


1.29 


1.80 


1.69 


1.76 


2.29 


1.71 


2.94 


1.65 


3.92 


1.53 


120136 


0.97 


0.18 


1.37 


0.17 


1.97 


0.17 


2.62 


0.16 


3.60 


0.15 


143761 


0.81 


0.89 


1.21 


0.87 


1.81 


0.85 


2.46 


0.82 


3.44 


0.76 


177830 


2.56 


4.19 


2.96 


4.10 


3.56 


4.02 


4.21 


3.90 


5.19 


3.65 


186427 


0.45 


6.59 


0.85 


6.44 


1.45 


6.29 


2.10 


6.08 


3.08 


5.65 


190228 


2.07 


8.11 


2.47 


7.91 


3.07 


7.71 


3.72 


7.42 


4.70 


6.82 


190360 


0.41 


18.81 


0.81 


18.36 


1.41 


17.93 


2.06 


17.33 


3.04 


16.06 


195019 


1.06 


0.53 


1.46 


0.52 


2.06 


0.50 


2.71 


0.49 


3.69 


0.45 


196885 


0.99 


4.23 


1.39 


4.15 


1.99 


4.07 


2.64 


3.95 


3.62 


3.71 


217014 


0.53 


0.20 


0.93 


0.20 


1.53 


0.19 


2.18 


0.18 


3.16 


0.17 



Note. — Values for were obtained from lCoxl <|2000h : Mx(G5V) = 3-5, M^(kov) = 3.9, Mx(K5V) = 4.5, Mk(mov) = 5.2, 
and M K(M5V) = 6.1. 
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Table 4. Calculated Secondary Star Angular Diameters 



Host Star Angular Diameter (mas) 

HD G5 V KO V K5 V MO V M5 V 



3651 





.77 


0. 


.71 


0. 


60 


0. 


.50 


0. 


23 


9826 





.64 


0. 


.59 


0. 


50 


0. 


.41 


0. 


19 


11964 


0. 


.25 


0. 


.23 


0. 


20 





.16 


0. 


07 


16141 





.24 


0. 


.22 


0. 


19 


0. 


.16 


0. 


07 


19994 





.38 


0. 


.35 


0. 


30 


0. 


.25 


0. 


11 


20367 


0. 


.32 


0. 


.29 


0. 


25 





.21 


0. 


09 


23596 





.16 


0. 


.15 


0. 


13 


0. 


.11 


0. 


05 


38529 





.20 


0. 


.19 


0. 


16 


0. 


.13 


0. 


06 


59686 


0. 


.09 





.09 


0. 


07 





.06 


0. 


03 


75732 





.68 


0. 


.63 


0. 


53 


0. 


.45 


0. 


20 


104985 





.08 


0. 


.08 


0. 


07 


0. 


.05 


0. 


02 


117176 


0. 


.17 


0. 


.44 


0. 


37 





.31 


0. 


11 


120136 


0. 


.55 


0. 


.51 


0. 


13 





.36 


0. 


16 


143761 





.49 


0. 


.45 


0. 


38 


0. 


.32 


0. 


14 


177830 


0. 


.15 





.13 


0. 


11 





.09 


0. 


.04 


186427 


0. 


.40 


0. 


.37 


0. 


31 


0. 


.26 


0. 


12 


190228 





.14 


0. 


.13 


0. 


11 


0. 


.09 


0. 


04 


190360 





.54 


0. 


.50 


0. 


42 


0. 


.35 


0. 


16 


195019 


0. 


.23 





.21 


0. 


18 


0. 


.15 


0. 


07 


196885 





.26 


0. 


.24 


0. 


20 


0. 


.17 


0. 


08 


217014 





.56 


0. 


.51 


0. 


44 


0. 


.36 


0. 


16 



Note. — The radii used in these calculations: G5 V = 0.92 Rq; 
K0 V = 0.85 Rq; K5 V = 0.72 Rq; MOV = 0.60 Rq; 
M5 V = 0.27 Rq. These were obtained from lCoxl l|2000h . 
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Table 5. Companion Check: Comparing Visibility Residuals 



AV max of Secondary Spectral Type Non Detection 



HD 


Obs Date 


0res 


G5 V 


K0 V 


K5 V 


M0 V 


M5 V 


Threshold 


3651 


2005/10/24 


0.043 


0.429 


0.481 


0.314 


0.187 


0.079 


M5 V 


9826 


2004/01/14 


0.024 


0.025 


0.018 


0.013 


0.007 


0.004 


G5 V 




2004/01/15 


0.022 


0.030 


0.021 


0.015 


0.008 


0.004 


G5 V 




2005/08/04 


0.102 


0.116 


0.095 


0.068 


0.047 


0.029 


G5 V 




2005/08/08 


0.063 


0.084 


0.064 


0.040 


0.022 


0.009 


G5 V 




2005/08/10 


0.065 


0.107 


0.086 


0.058 


0.036 


0.016 


G5 V 




2005/08/14 


0.060 


0.084 


0.064 


0.041 


0.022 


0.009 


G5 V 




2005/08/18 


0.051 


0.089 


0.068 


0.044 


0.024 


0.010 


G5 V 




2005/08/19 


0.055 


0.150 


0.137 


0.123 


0.125 


0.137 


- 




2007/09/05 


0.061 


0.031 


0.022 


0.015 


0.008 


0.004 


G5 V 




2007/09/06 


0.062 


0.034 


0.024 


0.016 


0.005 


0.005 


G5 V 




2007/09/12 


0.052 


0.042 


0.030 


0.020 


0.011 


0.005 


K0 V 


11964 


2005/12/16 


0.033 


0.204 


0.145 


0.086 


0.050 


0.020 


M0 V 




2006/10/20 


0.070 


0.222 


0.160 


0.095 


0.054 


0.022 


K5 V 


16141 


2005/12/12 


0.250 


0.334 


0.243 


0.148 


0.083 


0.033 


G5 V 


19994 


2005/10/27 


0.034 


0.229 


0.164 


0.098 


0.055 


0.022 


M0 V 




2005/12/10 


0.037 


0.229 


0.164 


0.098 


0.055 


0.022 


M0 V 


20367 


2005/12/12 


0.042 


0.402 


0.300 


0.190 


0.114 


0.050 


M5 V 


23596 


2007/09/11 


0.050 


0.278 


0.207 


0.127 


0.071 


0.030 


M0 V 




2007/09/14 


0.032 


0.270 


0.201 


0.123 


0.069 


0.029 


M5 V 


38529 


2005/12/06 


0.055 


0.081 


0.057 


0.033 


0.018 


0.007 


G5 V 




2005/12/14 


0.155 


0.066 


0.046 


0.027 


0.015 


0.006 


G5 V 


59686 


2005/12/06 


0.007 


0.007 


0.005 


0.003 


0.002 


0.001 


G5 V 




2005/12/16 


0.038 


0.008 


0.005 


0.003 


0.002 


0.001 


G5 V 




2007/04/02 


0.011 


0.007 


0.005 


0.003 


0.002 


0.001 


G5 V 


75732 


2007/03/26 


0.018 


0.101 


0.109 


0.103 


0.077 


0.037 






2006/03/30 


0.053 


0.201 


0.187 


0.149 


0.101 


0.046 


M0 V 


104985 


2006/05/17 


0.027 


0.009 


0.006 


0.004 


0.002 


0.001 


G5 V 




2007/04/26 


0.028 


0.009 


0.006 


0.004 


0.002 


0.001 


G5 V 


117176 


2006/05/13 


0.052 


0.277 


0.200 


0.120 


0.067 


0.026 


MO V 




2006/05/20 


0.019 


0.277 


0.203 


0.124 


0.071 


0.028 


MO V 




2007/04/02 


0.094 


0.277 


0.203 


0.124 


0.070 


0.028 


K5 V 


120136 


2005/05/12 


0.071 


0.018 


0.017 


0.011 


0.008 


0.005 


G5 V 




2006/05/14 


0.026 


0.040 


0.033 


0.021 


0.014 


0.007 


G5 V 




2007/02/05 


0.020 


0.001 


0.004 


0.004 


0.004 


0.003 


G5 V 




2007/03/26 


0.071 


0.002 


0.002 


0.002 


0.003 


0.003 


G5 V 




2007/03/30 


0.034 


0.001 


0.005 


0.004 


0.004 


0.003 


G5 V 


143761 


2005/06/29 


0.049 


0.250 


0.178 


0.107 


0.061 


0.026 


Mo V 




2005/07/03 


0.078 


0.251 


0.176 


0.103 


0.057 


0.023 


K5 V 




2006/05/12 


0.035 


0.251 


0.177 


0.105 


0.059 


0.025 


MO V 




2006/05/12 


0.092 


0.250 


0.176 


0.103 


0.056 


0.023 


KO V 


177830 


2006/08/08 


0.040 


0.101 


0.071 


0.041 


0.023 


0.009 


KO V 



Table 5 — Continued 



AVmax of Secondary Spectral Type Non Detection 



HD 


Obs Date 


<?"res 


G5 V 


K0 V 


K5 V 


M0 V 


M5 V 


Threshold 




2006/08/13 


0.095 


0.096 


0.068 


0.040 


0.022 


0.008 


G5 V 


186427 


2006/08/13 


0.072 


0.478 


0.358 


0.225 


0.130 


0.058 


M0 V 




2007/09/12 


0.048 


0.496 


0.371 


0.233 


0.134 


0.059 


M5 V 


190228 


2005/07/01 


0.076 


0.148 


0.106 


0.062 


0.034 


0.014 


G5 V 




2005/08/19 


0.044 


0.159 


0.113 


0.067 


0.037 


0.015 


K5 V 




2006/08/14 


0.042 


0.141 


0.102 


0.060 


0.032 


0.014 


K5 V 


190360 


2005/08/11 


0.045 


0.560 


0.421 


0.266 


0.156 


0.067 


M5 V 




2006/08/11 


0.044 


0.517 


0.394 


0.251 


0.147 


0.063 


M5 V 


195019 


2005/08/11 


0.049 


0.192 


0.149 


0.093 


0.054 


0.021 


K5 V 




2005/10/23 


0.086 


0.284 


0.212 


0.130 


0.075 


0.030 


K5 V 




2006/08/06 


0.070 


0.255 


0.191 


0.117 


0.068 


0.027 


K5 V 


196885 


2005/10/29 


0.053 


0.366 


0.269 


0.164 


0.094 


0.039 


M0 V 




2006/08/07 


0.037 


0.365 


0.268 


0.164 


0.093 


0.039 


M5 V 




2006/08/14 


0.055 


0.340 


0.252 


0.155 


0.088 


0.036 


M0 V 


217014 


2005/08/12 


0.047 


0.126 


0.100 


0.066 


0.040 


0.018 


K5 V 




2006/08/12 


0.033 


0.012 


0.010 


0.004 


0.0005 


0.002 


G5 V 
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50 100 150 200 250 300 550 

BASELINE (m) 

Fig. 1. — Example of the difference between the visibility curves for a single star and a binary 
system. The dotted line indicates the curve for a single star with 0=1.0 mas, while the solid 
line represents the curve for a binary system with the following parameters: pr imary=l-O mas, 
#secondary=0.5 mas, a=10 mas, and AK=2.0. 
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Fig. 2. — Example of residuals indicative of a single star. The solid line represents the 
theoretical visibility curve for a star with the best fit diameter, the dashed lines are the la 
error limits of the diameter fit, the Os are the calibrated visibilities (V c ), the vertical lines 
are the errors in V c , and the Os are the residuals to the diameter fit. Note how the residuals 
fall evenly both above and below the line at V = and show no sinusoidal trends. (Data for 
HD 120136.) 
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Fig. 3. — Example of residual systematics indicative of a binary star. The solid line represents 
the theoretical visibility curve for a star with the best fit diameter, the dashed lines are the 
la error limits of the diameter fit, the Ds are the calibrated visibilities (V c ), the vertical lines 
are the errors in V c , and the Os are the residuals to the diameter fit. Note how the residuals 
are systematically low for the shorter-baseline observations and systematically high for the 
longer-baseline observations. (Data for (3 Aurigae.) 
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2 4 6 8 

HOST STAR ABSOLUTE V MAGNITUDE 

Fig. 4. — Companion check: Eliminating companion spectral types. The x-axis represents 
the absolute magnitude for the observed exoplanet host stars and the y-axis represents the 
lowest-mass secondary that is still a possibility for each observed star. 
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Fig. 5. — Example SFP for a AK=0 binary system: (a) shows the primary star's fringe, 
shows the secondary star's fringe, and (c) shows the combination of the two. 
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Fig. 6. — Example SFP for a Ai^=2.5 binary system: (a) shows the primary star's fringe, 
(b) shows the secondary star's fringe, and (c) shows the combination of the two. 
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Fig. 7. — Example fringe envelope for a single star. The data for this star shows no indication 
of a secondary peak in the fringe envelope. (Data for HD 16141.) 
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Fig. 8. — Example fringe envelope for a separated fringe packet binary. The secondary peak 
represents the stellar companion. (Data for HD 3196. Image courtesy of Deepak Raghavan.) 



